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Abstract 


The definition of the geoid in view of different height systems is discussed. 
A definition is suggested which makes it possible to take the influence of the un- 
known corrections to the various height systems on the solution of Stokes' problem 
into account. 

A solution of Stokes' problem with an accuracy of 10 cm is derived which 
allows the inclusion of the results of satellite geodesy in an easy way. In addition, 
equations are developed that may be used to determine spherical harmonics using 
altimeter measurements, considering the influence of the ellipticity of the refer- 
ence surface. 
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1. Ihtrodttotion 


As part of its Earth and Ocean Physics Application Program (EOPAP) 
the National Aeronautics and Space Administration (NASA) plans in the next ten 
years the launch of some satellites equipped with altimeter for ranging to the 
ocean surface. The announced accuracy of the future altimeter systems lies in 
the scope of 10 cm. 

By solving the inverse problem of Stokes, it is possible to compute 
gravity anomalies from these very accurate altimeter data. In view of an ejami- 
nation about possibilities, problems and accuracy of a solution of the inverse 
Stokes' problem with this high accuracy we will treat two preparatory problems 
concerned with the direct solution. 

In order to transform the altimeter data into geoidal undulations' (if pos- 
sible, taking oceanographic informations about the so-called sea surface topo- 
graphy into account) we need d, suited definition of the geoid at least with the same 
accuracy. Another problem is concerned with the impossibility of the measure- 
ment of reasonable altimeter data on the continents. So we have to cut the dis- 
tant zones off in our integral solutions, taking into account their influences by a 
set of harmonic functions. This is also in agreement with the recent numerical 
treatment of Stokes formula (Vincent and Marsh, 1973, Bapp, 1973). The best 
set of harmonic coefficients is, of course, a combination solution. So we should 
have regard to the fact, that these coefficients do not belong to the potential on 
the ellipsoid or the earth surfece but to a sphere. 

The first comprehensive study of the direct solution of Stokes' problem 
with regard to the use of altimeter data is due to Mather (Mather 1973, 1974). 
However, because of the necessity of the inclusion of satellite coefficients into 
the solution, we will follow another way which seems better suite’d in'the. case of . 
our preconditions. 


The following treacmem: is oaetju on tntj robtutH uoacriucu m ^viurica, 
1974). In order to avoid long-winded repetitions we refer often to results and 
formulae of this study, so that the knowledge of this report may be recommend- 
ed for an entire insight into the present work. 

Eegarding our two special problems described above we first try to 
give a suited definition of the geoid. Of course, we shall not change the defini- 
tion of the geoid as an equipotential surface of the earth. What we will do is 
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nothing else than a specialisation of one equipotential surface distinguished as 
the geoid. Our main condition in this context is the possibility of a realisation 
of this special equipotential surface by geodetic measurements. 

A certain modification of Moritz's approach seems to be necessary if 
we want to include satellite data (e.g, in form of a set of harmonic coefficients). 
This is not so in view of the correction terms to Stokes' formula. But the in- 
clusion of satellite information directly in "Stokes approximation" (Moritz 1974, 
f. (1-9) ), may lead to very .complicated problems. 

In order of a better understanding of our problems and also the way 
which is chosen for the solution we will remember some basic considerations of 
geodesy. The main task of geodesy is the estimation of the figure of the earth 
and the outer gravity field with the aid of suited measurements. In the over- 
whelming cases these measurements belong to the earth's surface. 

We will assume that the earth's surface is a star-shaped surface. In 
this case any ray froin' the origin (the gravity center of the- earth) intersects 
this surface only once. Describing the physical surface of the earth by geodetic 
coordinates, the ellipsoidal height, can be considered as a function of the two 
other coordinates 


h = h(<P, X) 


h. ...... . ellipsoidal height 

geodetic latitude 

X geodetic longitude 


We will assume further that we have measurements of the following type 

a) ellipsoidal heights h (e. g. by altimetry over the- - 
ocean surfaces) 

b) potential differences C = W - Wo (by levelling) 

c) gravity g 
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Now by a combination of various t 3 ^es of these data we can obviously 
solve our main task in different ways. Because of the superficial similarity 
with the well-known boundary value problems of potential -theory we can also con- 
sider three geodetic boundary- value problems : 

a) First geodetic boundary value problem: Given the ellip- 
soidal height h and the potential difference C. Today, 
this method looks somewhat artificial, but with the re- 
cent development of doppler measurement methods or 
perhaps with altimetry on the continents, it may become 
very interesting. 

b) Second geodetic boundary value problem: Given the ellip- 
soidal height li and gravity g. The existence and unique- 
ness of the solution is discussed by (Koch and Pope, 1972). 

c) Third geodetic boundary value problem: Given the poten- 
tial difference C and gravity g. This is the well-known 
Molodenskii problem. For a detailed discussion of a so- 
lution see (Meissl, 1971). 

Supporting on this classification we will now make a few general com- 
ments, including a summary of some results. 

m all three problems we need a potential value Wq as additional infor- 
mation. It may be pointed out that by the inclusion of one additlonahpiece of data, 
(that is in case three, the inclusion of one geometric distance e.g, one ellip- 
soidal height h) the value Wo can be computed. 

It is nowadays impossible to measure C on the ocean surface. So the 
determination of the sea surface topography with the aid of geodetic measure- 
ments can only be obtained by the solution of problem two (Moritz, 1974), In 
the ease of the inverse problem we must assume that the altimeter information 
can be corrected by oceanographic information for sea surfe.ee topography, ha 
this case the corrected altimeter data should belong to an equipotential surface. 
Because of this assumption the information C = 0 is given in addition to the 
ellipsoidal height. If the equipotential surface is identical with the geoid, the 
ellipsoidal heights are identical with the geoidal undulations. 

The problem of the unknown geoid and the estlaaaation of datum parame- 
ter to the various heiaht svstems can be solved by combining data of all three 
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types. This is what we are going to do in the next two sections. The geoid as 
the. equipotential surfe.ce with the' value Wq is defined in such a way, that the 
square sura of the differences to the main height systems (corrected by oceano- 
graphic information about sea sirrface topography) is a minimum. To solve the 
problem of the practical determination of Wg and to compute the datum correc- 
tions to the height systems, condition equations of a least square adjustment' pro- 
cedure are derived in section three. It may be pointed out, that the "geoidal un- 
dulations", which are needed in this model as measurements, are not the true 
geoidal imdulations but values obtained from gravity anomalies which are falsi- 
fied by an unknown correction to the height datums. 

The main task of modem geodesy is not the solution of one of the three 
boundary value problems but a uniform solution which combines data of all types. 
From a practical point of view the combination of "terrestrial data" (gravity 
anomalies, altimeter data) and "satellite data" (orbital analysis, satellite to 
satellite tracking, etc. ) is the most important problem. At least the lower 
harmonic coefficients will be computed from a combination of all these data. 
From this point of view it is uncomfortable to use solutions of a boundary pro- 
blem, because the surface of the earth is very complicated. So it is important 
that the analytical eentinuatlon of the petentiai Inside the earth Is passible with 
any wanted degree of accuracy (Krarup, 1969). Moreover, if we start fronri 
the same data set (i.e. gravity anomalies) a series evaluation leads to the same 
formulae as the Molodenskii series solution, as shown by Moritz (Moritz, 1971). 


A computational procedure in which we can include terres trial anc 
satellite data is the following successive reduction method. 

A) Direct effect of atmospheric gravity reduction . Remove 
the atmosphere outside the surface of the earth and redistribute it 
inside. The resulting disturbing potential is then an analytic func- 
tion outside of the earth!s surface and the reduced gravity anoma- 
lies Agg are boundary values at the earth’s surface, 

B) Direct effect of the regard of topography. Compute gravit 3 
anomalies Ag^ at the geoid (or the ellipsoid) by a suited form of ana' 
lytical downward continuation, 

C) Direct effect of elliptic ity correction. Compute from the 
gravity anomalies: Ag; at the ellipsoid gravity anomalies Ag, at the 
sphere with the radius a. 
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At this state we can combine the gravity anomalies Ag; with satellite 
derived data. From the combined data we can compute the disturbing potential 
Ta at points on this sphere. We may remember that this potential must not 
be the true value of the disturbing potential at this point in space but only the 
result of the analytical continuation (consider the case of a mountain at the 
equator). 


D) Indirect effect of ellipticity correction . Compute the 
potential Te at the ellipsoid from the potential values T* at the 
sphere with radius a. 

E) Indirect effect of the regard of topography. Compute 
from the potential values Te at the ellipsoid the potential Ts at 
the earth’s surface by an upward continuation, using the inverse 
method of step B, 

i? ) Indirect effect of atmospheric gravity reduction. Cor- 
rect the value Ts at the earth surface by the indirect effect of the 
atmospheric reduction made in step A. 


Here, we will make only some remarks about this method and the re- 
sults. A detailed description together with a compilation of the formulae is 
given in section four. 

The .estimation of the direct and indirect effect of atmospheric gravi- 
ty reduction is the same as in (Moritz, 1974), The treatment of the influence 
of the topography is also very similar as used by Moritz. It can be shown ■ 
(Moritz, 1971) that the common handling of the direct and indirect effect leads 
to the same formulae as recommended in (Moritz, 1974, sec 4). However, 
the meaning of the procedure is quite different from Molodenskii's solution, 
which avoids analytical continuation. But the computational formulae are the 
same and well suited for practical computations. 

The treatment of the influence of the ellipsoidal shape of the reference 
surface is different from the procedure in (Moritz, 1974), The first part, the 
computation of gravity anomalies Ag, on the sphere with radius a from gravity 
anomalies Ag£ on the ellipsoid was done in the main in (Lelgoraann, 1972). 

The indirect effect, the computation of Te at the ellipsoid from T^ at the 
sphere, is derived in this report. 


- 5 - 



The use of the final formula for the correction term may also be fa~ 
vorable within the computation of spherical .harmonics from given altimeter 
data. Let us assume that we have altimeter data, corrected for sea surface 
topography, as a function of the geodetic coordinates. Then we obtain the dis- 
turbing potential at sea level by 


Te (<P, X) = N/y. 


N. geoidal undulation (from altimetrjt) 

y normal gravity 

Te disturbing potential at fhe'geoid or 

ellipsoid. 


In order to obtain a set of spherical harmonics we need the disturbing 
potential at the sphere with radivis a. 



We get this value by the simple correction 


($, X) = . Tg (9, X) - oos^ 9 • T {(p, X) 
e'p second eccentricity . 

The derivation of this simple formula is done in an indirect way. Be- 
cause of the length of the derivation it is given at the end of the report in the 
sections 5 to 6. 
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2. Considerations on the Definition of the Geoid 

The first definition of a geoid as that real eqmpotentiai suriace oi tne 
earth gravity potential, which is characterized by the ideal surface of the oceans, 
was given by Gauss. Such a definition presupposed that the ideal surface of the 
oceans is part of an equipotential surface of the earth gravitational field. To a 
■ certain degree of approximation this idealized sea surface coincides with another 
more or less time invariant conception^ the mean sealevel. 

We will consider here as mean sea level the mean ocean surface after 
removing time dependent effects. Because the mean sea level is than not neces- 
sarily an equipotential surface of the earth gTavitational field, slopes of mean sea 
level were detected both by' levelling and by oceanographic computations. The de- 
finition of an ideal surface of the oceans and the computation of the, difference be- 
tween this ideal surface and mean sea level caimot i)o a problem of geodetic but 
of oceanographic science. 

From a geodetic point of view the idea of a geoid is closely concerned 
with the definition of the heights. It is well-known that the heights are computed 
from measured potential differences. Let us assume for the moment the. (of 
course unrealistic) possibility, that we can carry oiit spirit levelling also over the 
ocean surfaces. In this case the geodetic community would certainly define as a 
geoid that equipotential surface, from" which the potential differences are counted. 

In every case the geoid must be considered as the reference surface of 
a world wide height system. So within the problem of the definition of a geoid 
the problem of the definition and also the practical possibilities of the computa- 
tion of height datums to the various height systems play a central role. 

There is a third utilization of the geoid or in this case rather the .quasi- 
geoid as established by Molodenskii, (Molodenskii, et, al 1962), which Is very im- 
portant in geodesy and this is the role of the geoidal undulations in the interplay 
of gravimetric and geometric geodesy. Apart from its own importance, we will 
use this connection to overcome the problem of the Impossibility of spirit level- 
ling over the ocean sur feces. 

The most important geodetic aspect in the considerations about a defini- 
tion of the geoid seems to be the definition of a reference surface for the height 
determination. So we will mention some principles which should be Important , 
in respect of our opinion that the computation of datum corrections is one of our 
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main problems. As in the ease of the definition of other coordinate systems 
(e, g. the definition of a highly accurate cartesian coordinate system for the pur- 
pose of the description of time dependent coordinates of geophysical stations) 
we may answer the following questions: 


1) What physical meaning has tne aennition r 

2) Can we transform the physical definition into a 
mathematical description ? 

3) Can we realize the mathematical and physical 
definition in the real world by measurements ? 

4) Can we compute transformation parameters to 
already existing height systems ? 


We will, give the ai^wers in the course of this section with exception of question 
three, which will be answered in the next section. 

After this preliminary considerations we will start the discussion with 
the definition of a geoid given by (Rapp, 1974). He started the discussion from 
the set of all equipotential surfaces of the actual gravity field 


W = W(x,y, z) = const. 


W is defined as the sum of the gravity potential Wg and the potential of the atmos- 
phere Wa. We must point out that the potential 


(2-1) 


•W = Wg + Wa 


is not harmonic -outside of the earth’s surface because of the presence of the 
atmosphere. 

. Because we are going to distinguish in the following considerations 


- 8 - 



several different equipotential surfaces, we will call an equipotential surface 
with the potential Wi {where the subscript i described only the fact, that Wj 
has a fixed value) 


geop (Wi). 


Later on we will specialize one of this equipotential surfaces as the geoid, that 
is 


geoid = geop (Wg). 


Departing from the customary ejq)ression for the potential of the geoid by Wo, 
we have characterized the potential of the geoid by Wg. The reason for the 
change of this abbreviation will be clearer in the course of this section. The 
choice of this special equipotential surface seems in a certain way arbitrarily, 
dependent on the starting point of the considerations. Por this reason we shall 
discuss, ,for the moment the problem separately from the three special areas 
we have mentioned at the beginning of this, section. Then we will look for a com 
bination of all these considerations. Most important, of course, is the possi- 
bility of a practical realisation of the geoid in the case when we have enough 
accurate measuring data. 


1) The geop (Whsi.) as the ideal surface of the oceans. 

(MSL mean sea level) 

The definition of mean sea level and the ideal surface of the oceans, 
which we will consider as an equipotential surface, cannot be the task of geodesy 
but of oceanography. A good description of the difficulties of the definition and 
more over the realisation of these concepts are given in (Wemelsf©lder.’, 1970), 

The following very simple model of the real processes may be sufficient 
for our considerations. Because geodesy is only interested in the deviation of 
the ocean surface from an equipotential surface, we may say, that the ideal sur- 
face of the ocean is disturbed by the following irregularities: 
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a) very short periodic irregdlarities (e. g, oceai 
waves, swell) 

b) periodic or quasi periodic irregularities (e. g, tides) 

c) quasistationery irregularities, which retain their form 
but change, their place (e.g. gulf stream) 

d) quasistationery irregularities, which retain form 
and place 


If we correct the real ocean surface for all these irregularities, than the- 
result should be an equipotential surface and we will name it by 


geop (Wmsl) 


Corrections of the individual height datums to a world wide height system 
are then given by the correction (dWMsu) for quasistationery sea surface topo- 
graphy at the water gauges. 

Such a definition is not only important for oceanographers, but also for 
geodesists. If we can compute with the help of oceanographic information the de- 
viation of the sea surface from an equipotential surface, then we can also com- 
pute geoidal undulations from altimeter measurements. Especially if we want to 
recover gravity anomalies from altimeter data we have to use such information. 


2) The geop (Whso) as the basis of a worldwide height system 
(HSO. ............ height system zero order) 

In order of an explanation of a geop (Whso) let us start froni a reference 
surface of a particular height system geop (Whsi) (e. g. from the mean sea level 
1966.9 at Portland, Maine). We will assume errorless levellings to the reference 
points of (n-1) additional height systems. Consequently, we have n different 
height systems with the reference points on equipotential surfaces 


geop (W„3i) 
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Of course, sucli kind of levelling is impossible because various height 
systems lay on various continents. We shall bridge this difficulty using the 
connection between gravimetric and geometric geodesy. 

The most plausible reference surface of a worldwide height tsys tern 
is then the equipotential surface geop (Wh 3 q) for which the sum of the sqpiare 
deviations to the particular height systems is a minimum. 


( 2 - 2 ) 


n 



w^f =Min^ 


In this case all height systems have equal Influence. As a solution of the pro- 
blem we get easily 


(2-3) 


Whso = 


_1 

n 


n 


I 


Wh31 


1 = 1 


We can assume that a geop (Whso) defined in such a way lies very near the 
geop (Wmsl), because all height systems are based on mean sea levels at least 
in the reference points. The transformation parameters are given by the defi- 
nition equation. 


3) The geop (Wq> from the connection between gravimetric and geometric geodesy . 

(Wo = Uq Uo = normal potential on the surface of the normal ellipsoid) 

We will start the definition of a geop (Wg) from the normal potential 
based on a rotational ellipsoid. The surface of the ellipsoid should be an equipo- 
tential surface of the normal potential. It is well known that in this case the nor- 
mal potential on and outside of the ellipsoid and also the geometric form of the 
ellipsoid itself can be described by four parameters, e. g. 


k M mass of the earth 

CO rotational velocity 

Jg harmonic coefficient 

of order two 
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Uo. ....... potential on the ellip- 
soid surface 


or 


a. semi major axis 


We presuppose that we have very exact values- of the first- three terms 
(maybe from satellite geodesy). To a certain degree 'of accuracy the following 
relations hold (Heiskanen-Moritz, 1967) for the fourth term 


(2-4) 


Uo 


kM 

a 


(2-5) 


dUo 



- yo^a. 


• It is well known that we cannot measure the absolute value of the poten- 
tial. We get this value by an indirect method using e. g. the connection between 
gravity and a distance in formula (2-5). For a definition of the normal gravity 

field, it is important that only three physical constants are fixed values of the 
real earth. The fourth term is in certain limits arbitrarily. Let us describe 
the geometrical relationship in this case. We have the following situation. 
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hp ellipsoidal height 

Np geoidal undulation 

Hfo .orthometric height of the 

point P in a worldwide 
height system 


It is easily seen that in the case outlined above, the main equation, • 


(2-6) 


h = H + N, 


connecting gravimetric and geometric geodesy holds not in this form. We have 
two possibilities to correct the situation. We can refer the heights to the re- 
ference surfe.ce 


geop (Wo) 


or we can change the size of the normal ellipsoid by 


d Uq L T ( Wuso W o) 


The choice of the'‘kind"of the‘'correelion is our own pleasure. On the '•“ 
oilier hand if we''have fix,ed"one of both" values ''((liat'is Wh3o""OX’.Wo» ‘eilher'by a> 
mark on "the earth sun"face'or by a-giveu number); i:he diflierenee between them'" 
must be computed frona geodetic measurements. 


In this connection we will also consider the mathematical, description, 
of our problem. For this purpose we use the relation 




I II 


h= =„r..+ih,o., 


hn . ellipspldal„hei,gh,t,.r„ 

H prthometric, height 

N . . HU . . geoidal.undulation.ii- 

. ij . . .iqiianioTgeoidali.iundulatloniH 

IP'. . w. .Eormal-heigM. 
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For an accurate definition and ejiplanation of all these terms see 
e.g. (Heiskanen- Moritz, 1967). 


4) The geoid - geop (.Wg) as a result of the previous considerations 

Because of the close connection with a worldwide' height system our 
final definition of the geoid is based mainly on point two of the previous con- 
siderations, In addition we will take into accoimt information from oceano- 
graphic science. In order to bridge the continents and to get a numerical 
value of Wg we shall use also-the connection between gravimetric and geome 
trie geodesy. We study the geometrical relationship on the following figure. 



The geoid is defined by 

n 

(2-8) X y " 


Hpi orthometric (or normal) height 

in a particular height system i 

dHfisi reduction of the gauge (mean sea 
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level) mark to an ideal ocean 
surface because of sea surface 
topography, given from oceano- 
graphic science 

dHHsi deviation of the corrected basic 

level of the height system i from 
the geoid, defined by (2-8) and 
constant for the area of this parti- 
cular height system 


da correction of the semi-major axis 

for the term (Wo - Wg) 

N geoidal undulation 


In an explicit form we have for the definition of the geoid 


(2-9) 


Wg 





dWwsi). 


If the oceanographic corrections dWHsi are correct then all geop 
(Whsi + dWMsi) are identical or at least close together and the co rrec tions dWnsi 
are small, or zero. In this case the practical procedure developed in the next 
section may be regarded as an independent checking of the oceanographic infor- 
mation by geodetic methods. If we have no oceanographic information we can 
put, in this case, simply dWMa ~ 0, 
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3. On the Realisation of the Definition of hhe Geoid 

In the previous section the possibility of transferring the theoretical 
definition into physical reality and vice versa was one of our main requirements. 
This is certainly a question which can only be answered by statistical methods, 
that is by the development of a suited adjustment model. 

Our mathematical description of the problem is based on the connection 
between gravimetric and geometric geodesy. So we can start with the condition 
that the basic equation of gravimetric geodesy 


h - H+N = C+H* 


is fulfilled in a set of m points. Because the observations from levelling, from 
geometric and from gravimetric geodesy may be given in different systems we 
must include in our model transformation parameters as unknowns. lii this way 
we are lead to the model of a least square adjustment of condition equations with 
additional unlcnowns 


(3-1) 


j?iv+I?x+w = 0. 


The solution of such a system is well known (e.g. Gotthardt, 1968, p. 238 ff). 
We will discuss here the explicit form of the condition equations, presupposing 
that the following "observations "are given at m points Pj on the earth surface: 


hj ellipsoidal height, computed from 

rectangular coordinates as a result 
of satellite geodesy 

Hfj normal height in the i-th system- of 

n height systems 

quasi-geoidal undulation 
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In addition, we have for any of the n height systems a constant dH«st, repre- 
senting sea surface topography as computed in oceanographic science: 


dH^si gauge correction due to sea surface 

topography in the i-th height system 


In order to fulfill the condition equation 


(3-2) h - H* - C = 0 


we will first consider our "measurements” (hj, , Cj). 

As mentioned above these "observations" may belong to different sys- 
tems. The transformation parameters between these systems may not be known 
and have to be estimated in the course of the adjustment. This is true in any 
case for the height datums dEtia and the correction of the semi-major axis da, 

In this way it is possible to take also other systematic effects into ac- 
count, We will restrict ourselves to the unknown parameters described in' the 
following context. 


a) Ellipsoidal Height h 

We assume that the geometric reference ellipsoid is not in an 
absolute position, but connected with the center of mass by the vector {dsy, 
dyo, dzo). Then for absolute ellipsoidal heights h we get the following equa- 
tion (Heiskanen- Moritz, 1967, p. 207) 


(3-3) h = hj - cos75cosXdxo - cos^sinXdyo - sinpdzo - da, 

We have already included in this equation the unknown 


(3-4) 
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from which we can compute the numerical value of the potential at the geoid bj 


(3-5) 


Wg - Uo + dUo 


after the adjustment. 


b) Normal Height H* 

From figure two we get easily 


(3-6)' 


H* = + dHrtsi + dUHSj 


The n unknowns d Hhsi are the transformation parameters for all the height 
systems together with the constant values dBisj*. 


' c) Quasi-geoidal undulations C 

This case is less trivial, because the "measurements" Cj are 
connected with the unknowns dEwsi by Stokes' integral formula. With regard 
to our adjustment model we must derive a linear relationship between the "mea- 
surement" Cj and its true value C 


(3-7) 


c 





dHnSk) 


The influence of by the unknowns dH^sk is due to the fact that we can compute 
gravity anomalies Ag only with heights related to the geop (Whsu + dWnsk). For 
gravity anomalies Ago related to the geoid we have the expression 


(3-8) 


Age - g - Te + 0. 3086 [(’Hjt + dHwsfc) + dHnskJ 


or 
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(3-9) 


ZigG = Agfc + 0.3086 dHnsk 


with 


= (g " ■yt + 0. 3086 (Hjc + dKtisk) )• 


We have used in this formulae the normal gradient of gravity 


F = 0.3086h 


obtaining F in mgal if we introduce h in meters. The aecurady of the formulae 
above is sufficient for the computation of the coefficients a^j. 

The use of Agi^ is in agreement with our conception, that the first step 
to the estimation of the geoid is the inclusion of oceanographic information. The 
gravity anomalies Ag^ are related to the geop (Wwsk + dWHsk), which should be 
very closed to the geoid or, in an ideal case, should already coincide with the 
geoid. If we put the expression (3-9) into Stokes* formula, we get 


C = Cl + f'C 


where 

(3-10) 


fC 


R 

4ttG 


(O.3086 dH«sk)S(!^)da. 


In spherical approximation we can write 


0.3086 = 


d.y ^ 2G 

rih R 
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Putting this relation into (3-10), we obtain 


(3-11) fC = S(^)da. 

n 

The area in which heights related to the geop (Wnsk) are being used, may be ■ 
described by FHSk* 

Because dEnsk is constant over the area Fwsk, we get the loiiowing 
expression 


fC - ^ dlW • Js(#da 


or 

(3-13) a^j = J s(j/?j)dCT 

Fnsk 

where S(^j). is Stokes’ function referred to the point Pj in question. From the 
expression (3-13) we are able to compute numerical values for all the coeffi- 
cients atj, replacing the integral by a summation. 

In order to get a feeling about the magnitude of ait we will consider a 
simple example. We assume only one height system with a.feference surface 
different from the geoid. We assume further that the area Fhsi of this height 
system is a spherical cap of size around the point Pj under consideration. 
Starting with 

00 - '2u 

1 r 

= — i — S(0) sin0d0dQ! 

2tt J " 

0=0 O!=0 
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the integration over cl gives 


a-kj = f S(^ti)sin^d0. = 2j(^o) 

^ i )=0 

where the function J(i/)) is defined by (Heiskauen-Moritz, 1957, p. 119) 


JUl)) = ^ f 

« 

0 


S(i/j) sin;/)d0. 


From a table of . (e. g. Lambert and Darling, 1936) we draw the 'values 

J(i/)) = 0.5 at 0 = 27'" and 0 = 50° . 

Hence, raider this circumstances the coefficient of the unlcnown dHnsi is 

(1 ai.j) = 2, 

After this preliminary discussion about the connection of the "measurements’* 
and the unknown parameters, we will now derive the eiqilicit' form of the 
condition equations. We put (3-3), (3-6) and (3-7) Into (3-2) 

(hj + Vij) - (H=^j + Vgj) - (Cj + Vgj) - cos(pcosXdxo 

( 3 _ 14 ) - cos(psinXdyo - sinipdZo - da - dBisi - dEwsi 



dHnsk) - 0, 
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To this system of condition equatioxis, we have to add the equation 


(3-15) 



which is a consequence of the definition equation (2-8), This additional condi- 
tion equation connects only unknowns and not the "observations". However, 
there are no principal problems concerned with the solution if we add equation* 
( 3 - 15 ). at the end of the^equation system (Gotthardt, 1967). 

We will now discuss some main aspects related to the present and fu- 
ture accuracy of the data. We seperate the discussion into four parts in accor- 
dance with the different type of data which are needed in the adjustment. 


a) Determination of gauge corrections dlkisi: 

For purposes of physical oceanography. to interpret the results of a 
satellite altimeter, it is necessary to Icnow the absolute shape, of. a level surface 
near mean sea level. Also for geodetic purposes it is necessary that we know 
the deviation of mean sea level from an equipotential surface. We can com- 
pute this deviation combining altimeter measurements and gravity measure- 
ments (Moritz, 1974, sec. 5). However, gravity measurements over sea are 
very time-consuming. So from a practical point of view it could be very help- 
ful if we can determine gepidal undulations directly from altimeter measure- 
ments taking a small correction term from oceanographic science into account. 

To compare the results from oceanographic science with the geodetic ' 
results it is of course very important to determine the heights of several gauges 
in the same height system. Differences. between theory and measurements arc 
known (AGU, 1974). The first step in the realisation of a worldwide geoid 
should be an ejq)lanation of these differences, because this seems to be an indi- 
vidual problem which can perhaps be solved prior to the inclusion of geometric ■ 
results and geoidal computations. From a comparison of levelling and the 
results of oceanographic science at the coasts of the United States (AGU-,- 197,4) 
we can conclude that the present accuracy of the oceanographic, computations 
of sea surface topography is better than a meter. 
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b) Determination of station coordipa(:es and the ellipsoidal height h,i 

We expect accurate station coordinates from satellite geodesy, deri- 
ving ellipsoidal heights from the cartesian coordinates. The present accuracy 
lies in the order of 5-10 m (Mueller, 1974). However, with the high accuracy 
of new developed laser systems (some cm) and with special satellites, such as 
Lageos or special methods, such as lunar laser ranging we can expect a fast 
increase of the accuracy at least for the coordinates of some geophysical sta- 
tions. So it seems advantageous to use such geophysical stations also as the 
basis for the ac^ustment model developed above. It should be useful to have 
several well distributed stations in the area Fhsi of every main height system. 
Necessary are likewise the very accurate parameter of a normal ellipsoid, that 
is at least the three parameters { co, Jg, 1<M), computing the semi-major axis-' 
a within the definition of the geoid. 


c) Determination of normal heights' H* or orthometric heights H 

For long times levelling was one of the most accurate geodetic measux- 
ments with a standard error up to ±0. 1 mm per km distance. Today, neverthe- 
less, the increasing accuracy of distance measurements let tis expect a similar 
accuracy for coordinates and distances. To stay comparable in the accuracy a 
very careful examination of systematical errors in levelling is required. 

So far it is possible, a mutual connection of the geophysical stations 
and also the connection of these stations and the fundamental gauges of the height 
system by high precision levellings should be performed. 


d) Determination of geoidal heists N or quasi-geoidal undulations C 

This seems from the present accuracy considerations the most crucial ’ 
point in the method. On the other hand the definition of the geoid given above is 
connected with precise gravimetric geodesy, so that we can assume the pre- 
sence of data with the necessary accuracy. 

This is not the case for present gravity data which allows a computa- 
tion of geoidal undulations with an accuracy of better than 10 m. Methods in 
recent development like aero gradiometry and so on will provide perhaps a much 
better estimation. 

However, a highly accurate determination of the harmonic coefficients 
of lower order and a good and dense gravity material around the geophysical 
stations seems of high value in the solution of our problem. 
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4. On the Solution of Stokes’ Problem Including Satellite Data information 

We will nof; consider Smokes’ problenn as <-he compu<^a<;ion of a regularized 
geoid.but as the following problem: It may be possible to determine the disturbing 
potential on and oUtside of the earth,- using gravity data from the earth's surface 
in Stokes' integral formula and correct the result by some small terms. In this 
way we .can" interpret Molodenskii's solution, as a well suited and theoretical unob- 
jectionable solution of Stokes' problem. 

On the other hand, we have to go nowadays a Step mrtner. uur uttta, 
comprise not exclusively gravity measurements on the surface of the earth but al- 
so" a lot of other information about the gravity potential, one of the most important 
the lower harmonic coefficients from satellite geodesy. 

We do not know an exact solution of Stokes' problem. We wiir discuss on- 
ly approximate solutions, but the approximation error must be less than 10 cm in 
the geoidal undulations and the solution should be as simple as possible in view of 
practical computations. 

Because of the presence of the atmosphere, we have to solve not a La- 
place but a Poisson equation. We shall overcome this difficulty by ‘a suitable . 
gravity reduction. 

The majority of the data (the gravity anomalies) calls for a solution of 
the so-called Molodenskii problem. This is a very complicated type of a non- 
linear boundary value problem, Molodenskii already has based his so lution 
on Stokes' formula, because the result of Stokes' formula differs from the exact 
solution only by small correction terms. 

A unified treatment of Stokes' problem with the necessary accuracy of 
better than 10 cm for geoidal undulations was done by Moritz (Moritz, 1974), 

Such a solution'must be taken into account, 

1) the effect of the atmosphere 

2) the influence of topography, 

3) the ellipticity of the reference surface. 
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Neglecting terms of Mgher order of a Taylor series eijqjansion, Moritz 
treats all three effects independently of each other. He ends with the following 
procedure; 


1) Reduce Ag to the "Stokes' approximation" Ag° by, 

3 

(4-1) Agp = Ag - y Gi. 

1 = 1 


2) Apply Stokes ' integral to Ago 

(4-2) JJa^ S(l/))dCT. 

a 

3) Correct to obtain the actual value C by, 

3 

(4-3) c = + T Z,. 

1 = 5 . 


This solution of Stokes' problem is well suited if only gravity anomalies 
are at hand. However, the inclusion of data other than gravity anomalies may im- 
prove the solution very much. At present, the most important additional data are, 
without question, the results of satellite geodesy in the form of potential ooeffir 
cients. E seems very difficult to include this data in formula (4-2) in a convenient 
way. 


For this reason we have to change the model of Moritz. For a better unde 
standing of the basic idea, we will explain the difference in the case of the applica- 
tion of the ellipsoidal corrections. 

In (Moritz, 1974) and also in CLelgemann, 1970), we have established a 
one to one mapping of the reference ellipsoid on a sphere- with radius R by mapping 
a point ,P of geodetic (geographical) coordinates (<p, X) on the ellipsoid into a point 
P' of. spherical coordinates ($ = <p,X) on the sphere. 
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In this case the sphere is not attached to the ellipsoid, serving only as 
an auxiliary surface for the computations and the result belongs, of course,, to the 
point P on the ellipsoid (or rather to the earth’s surface). ■ 

Now, we will describe the present model. In this case, we compute from 
the anomalies on the ellipsoid the anomalies on the sphere with radius, a, which is 
tangent to the ellipsoid at the equator. With the aid of Stokes’ integral, we com- 
pute Ta at the point P’ on the sphere (if P’ lies outside of the earthlslsurface, 

Ta is the real disturbing potential at the point P' in space). At this step, we can 
combine the surface data with satellite data. Finally, we compute from the dis- 
turbing potential at the sphere with radius a the disturbing potential Te at the sur- 
face of the ellipsoid. This treatment of the ellipticity seems to be the appropriate 
ejq)ansion of the correction for spherical approximation in view of the fact, that the 
analytical continuation of the disturbing potential is possible with any wanted accur- 
acy. 


In order to get a closed theory, we must also treat the influence of topo- 
graphy in another way than Molodenskii. Molodenskii's solution is identical with 
the amlytical continuation to point level (Moritz, 1971). The formulae for ana»- 
lytical continuation are derived by Moritz in the cited publication. Because there 
is no theoretical difference in the reduction to different level surfaces (Moritz, 
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1971, sec, 10) we can first reduce the measurements to sea level and afterwards 
hack to the earth's surface. The formulae remain nearly the same, but we should 
look at the terms of &e present model as a result of a computation of the potential 
of analytical continuation at the ellipsoid and the sphere with radius a. In this 
way we end with the following procedure: 


1) Reduce the gravity anomalies for the effect of the atmospher' 
(4-4) Ag, == Ag - 6gi. 


(4-5) 


2) Compute by a suited form of analytical downward continuation 
gravity anomalies at the geoid or better, immediately at the 
ellipsoid 

AgE = Ags - 6ga. 


(4-, 6) 


3) Compute from the gravity anomalies on the ellipsoid the gravity 
anomalies on the sphere with radius a 

Ag^ = AgE - 6gs . 


(4-7.) 


4) Apply Stokes' integral in order to get the potential (of analytical 
continuation) at the sphere with radius a 




4n 


JJ (Aga) S(^/))dCT. 


a 


(4-8) 


5) Compute the potential at the ellipsoid by 


Te = Ta + etg . 


(4-9) 


6) Compute the disturbing potential at the earth's surface by up- 
ward continuation 

Ts = Te + 6t, 
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7) Correct the disturbing potential at the earth’s surfece by the 
indirect effect of removing the atmosphere 

(4-10) T = Ts + 6tT . 


From the disturbing potential T we get immediately the quasi geoidal 
undulations by 


(4-11) C = 


and tbe geoidal undiilations by 


(4-12) 


N = C + (H* - H). 


Now, we will collect the formulae for the correction .terms, together 
with the references. For a detailed explanation of these formulae see the re- 
ferences. 


(4-13) 1) figa. = - - = - 6gA . 

(See Moritz, 1974, formula (2-23) ). Note that this cor- 
rection has the same absolute value as the gravity correc- 
tion 6gA in (lAG, 1971, page 72). 1 takes its maximum . 
value of 1 6gA- I = 0. 87 mgal at sea level. 


(4-14) 2) 6ga =; H • Li (Ag) 


with 

(4-15) 


I^(Ag) = 




2tt 


r -M. da. 


- 28 - 



(4-16) 


(4-17) 


(4-18) 


(4-19) 


(See Moritz, 1971, formxdae (1-5), (1-8) ). Note th? 
this is only the first term of a series solution. . How- 
ever, it seems to be sufficient for all practical purposes, 
especially in the case of the use of altimeter data). 

3) ^ ^ (n-1) [C„«R„«(Q,X) + D„„S^„(e,X)] 

n =St>n — 0 

The coefficients can be computed in the following manner. 

When 


00 ^ 

T(0,X) = ^ ^ [An^„«(e,X) + Bu»S„„(0,X)l 

n=ara=0 


then 


Cna + Aj,n Qan + 


Dtto BnaQras B^Q)o^nm 


where. 


Qffl 


Qna 


(5n - 17) (n - m - 1) (n - m) 

4(n‘- I) (2n -■3).(2n - 1) 

-6n^ + 8n^ + 25n + 6nm^ + 6m^ + 21 
4(n - 1) (2n + 3) (2n - 1) 


• no 


(5 n + 11) (n -H m + 2) (n -f m + 1) 
4(n - 1) (2n + 5) (2n + 3) 


See (6-33), (6-35) and (6-41) and also (Lelgemann, 1972, 
formula (33) ), The formulae give the reduction term in 
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the case of gravity anomalies. Similar formulae for 
the computation of gravity disturba.nces at the sphere 
from gravity disturbances at file ellipsoid are derived 
in section 6 of this study. 

JI S(i/))dor. 

a ■ 

The disturbing potential T is computed instead of ge- 
oidal undulations in order to avoid the definition of 
geoidal undulations or quasi geoidal undulations in 
space. 


(4-20) 


5) 



T(6.X). 


The derivation of this simple formula is rather lengthy. 
It is given in the last three sections of this study, (see 
formula (6-35) ). 


(4-21) 


(4-22) 


(4-23) 


6) 6{^ = ■- H • Ag 

(See Moritz, 1971, formula (1-14) ). Note that this _ 
correction is almost zero in the case of altimeter data 
because of H = 0. 



r 


or 

'co 

fiti = J 6gi(r)dr’ 
r 

(See Mpritz, 1974, formulae (2-22) and (2-27). Note 
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that the indirect effect amounts to maximal 0, 6 cm. . 
hi spite of the overall accuracy of the solution, we 
can neglect its influence. 


All the correction formulae are the same in the case of gravity dis- 
turbances with the exception of the correction term Sgg.- The correction term 
6g3 for gravity disturbances is given by formula ( 6 - 33 ). 
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5, The Influence of the Ellipticity in the Case of Gravity Disturbances, 

The purpose of the next three sections is a unified treatment of the in- 
fluence of the- ellipsoidal reference surface and an unified evaluation of the formu' 
lae which connect the potential on the ellipsoid and the potential in space (that is 
on the sphere with radius a). The content of the present section shoiild only be 
seen as an intermediate result, which will be needed in the following section. 

In contrast to the very similar derivations in (Lelgemann, 1970) the 
derivation of the whole theory is based on gravity disturbances. The advantages 
rest on the avoiding of the difficulties concerned with the spherical harmonics of 
zero and the first order, which appeared in (Lelgemann, 1970). 

We will mention that the first explicit solution based on gravity distur- 
bances as data was derived by Moritz (Moritz, 1974), starting from the solution 
of the problem for gravity anomalies. Using his technique in an inverse way we 
shall derive in the next section the solution for gravity anomalies from the solu- 
tion based on gravity disturbances. 

On the surface of the normal ellipsoid we have in a linear approxima-r 
tion the following boundary condition. 


(5-1) 



= -6g. 


Using Green's second formula for the function T and the ellipsoid as the integra' 
tion surface. 


(5-2) 


-2ttT 





dE 


we get, after inserting the boundary condition, the integral equation 
(5-3) . jT(-^(-^))dE. 
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In consequence of the assumption above, this is a Fredholm integral equation 
of the second kind. The integrals must be taken over the ellipsoid. In order 
to get a solution we transform first the integrals into a spherical coordinate 
system. The result is a linear integral equation with an unsymmetric kernel. 
We shall develop this kernel with the required accuracy into a power series of 
e'®. The resulting system of integral equations consists only of equations with 
symmetric kernels, moreover of equations with a well known kernel.- Because 
the eigen-functions of the integral equations are the spherical harmonics we 
arrive very easy to series solutions. 

We map points of the dllipsoid with the geodetic latitude in such a 
manner onto the sphere that the spherical latitude is identical with the geo- 
detic latitude (left side of the figure) 



0 = 90° - <p 

First we shall evaluate some terras in the powers of e’® up. to the re- 
quired degree of accuracy. For the surface element of the ellipsoid we have, 

(5-4) M = a(l+e’®)"^^ (l + e'®cos®(p)“ '^= a(l-e'® 3/2e'®cos®6) 

(5-5) N = a(l-fe'®)^/® (l+e'®cos®<p)“'*’'^ = a(l -tie'®cos®0) 

(5-6) dE = M-N* cos<pd<pdX = M*N*da = (l-e'®+ 2e'®cos®6) a®dq. 
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Further we will represent 1/4 in dependence of- 1-/^ (see figure above). -We 
start with 


(5-7) 4® = (3Q- - -52)^ + (Ye - Ye')^ + (Ze - 


where, 

Xe = Ncos¥>cobX, Ye = Ncos^osinX, Ze = b^a^Nsin^O. 
Developing this -formula in powers of e’® and using the substitution. 

(5-8) = (X - X’)® -f (Y - Y’)® + (Z - Z'f 

with 

X = acospcosX, Y = acosv5sinX, Z = a sin<D 


we get 


(5-9) 4® = + -|e'^(co3^0 - eos®9’) - 2e’?a^(cos6- cos 9’)' 


and with the help of a series evaluation the final result, 

(5-10) - ~~ e’®(cos®9 + cos®9') + e’^a® ^ ' J 


An exact expression for — 2 — (1/4) was derived by Molodenskii (Molodenskli, 
1961, page 54) 
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( 5 - 11 ) 


= 


— (1 + e’®n®) 


( 5 - 12 ) 


1 

-n = — 5 " (Nsin<p - N’sin<p’). 

JC£ Si 


A series evaluation gives 


( 5 - 13 ) 


^ m i - _1_ fi . 4 e'^oos^e H. 6'=' (0°s9 - cos6') l\ 

5n \ 4 / 2a4 V • ■' 


( 5 - 14 ) 


= 1 /a • i = 2 sin (&/ 2 . 


Inserting these expressions into the' integral eq\iatieh-( 5 - 3 ) , we obtain, after 
some transformations 

2 nT = J -[l-e'® +-^ e'^cos^e- -~'e'®cos®e' + 

( 5 - 15 ) • fig • _S_da + r - _^e’®cos *0 + _i_e'®cos^ 6 * 

^ L 2 2 8 8 

- e'^ .(cosQ- cose’)^ -I . T . J-. da. 

Iq J ■^o 

The kernel of this integral equation is not symmetric. So we evaluate T in 
powers, of the small term e'^ 


( 5 - 16 ) 


T = + e'^T^ + e'^T^ 
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Inserting, this series into the integral equation (5-15) and equating the coeffi- 
cients of the powers of e'“ we arrived at the following integral equation system 


2ttT = 


= [a • 6g • — dcr - -g — 

I ^ ^0 I ^ 


dCT 


2ttT'" = 


* =|a.. 6g. + ^ 


-L. cos^e - -i- cos*e» 

4 


(5-17) 


da + 

r 




COST'S 


1 ^ - ooBWf ] aa - ^ da. 

8 Ir, J 2 J 


iTlX'' = 


we musi. solve successively the first and second equation of this’ system. These 
two equations contain all terms up to the order e*® so that a solution of the fur- 
ther equations not necessary in view of the required accuracy. 

We develop the first integral equation into a series of eigenfunctions. 
Because, 


(5-l») J X„(0,A) P„(cos0)da = Xn(0',X«) 

O' 


we get the series representation. 


(5-19) 


T 



JSs_ 

n + 1 
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Analogous to Stokes formula we can represent this series as an integral formu- 
la (Hotine, 1969) 


(5-20) 


T 


a 

4tt 


J fig S(0)dcr 

CT 


With 


(5-21) 


^ f P.(COB0) 

Li n + 1 

n = 0 


or by the closed expression 


(5-22) 


S(0) = 


sin 0/2 


- f 1 + 


sin 0/2 


')• 


In order to get a solution of the second equation of the system (5-17), we evalu- 
ate the disturbing term in a series of harmonic functions. Because is to be 
multiplied with e’** we can use the following formulae in spherical approxima- 
tion: 


(5-23) fig = -i- y (n + 1) T„ 

^ tl= 0 


and 


(5-24) 



CO n 

{Ann Rna(®» Ban ®nn(®» 

n = u n = o 



As an intermediate result we get 
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{5-2E 

" 4. X. ° " 

As I Afl + A.j 

Wltll 



(5-26) 

- L ] 

or n:=Q 

= -Vt„ 

■^0 J^O' 

(5-27) 

O’ ” 

cos'® 6 ' V <71 
4 " 

Using formula (8-6) we obtain 



. ^ 1 ’ r (cos6 - cos6*)^ ^ , 

- ir, 1 4 

2* n* Tft dCT = 

(5-28) 

f f 2n(4m^- 1) ■ 

L L (2n + 3)(2n-l)(2nH 1) 

n=0iB = 0 


and with formxila (8-5) 




1 2 q da 

cos 0 

J • S,Q 

(5-29) 

n 

4 r («■'"«> I [^“{^ 

ISO B =0 

» (0.^) 




2n + l 


Roo (®» X) + 


2n-3 


R ^-a) 
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l [>^8 



+ B 




{ 




2n + 5 


S ^+2) n(®» 


s, 




2n + l 


Sum (®» “*■ 


y, 


iUl- 


2n-3 


'r») 


..(e.X) }]. 


Performing the integrations on the left hand side of the integral equation we have 


(5-30) 


4. 1 r da ^ ^ n + 1 rp i 

• 4tt J JIq 2j 2n + l ” * 

a n =0 


Transforming first the term Ag in a series of harmonics we get afterwards from 
a series comparison on both sides of the integral equation the following solution, 


09 


2n + l 
n + 1 


,U 


nfW - 1) 

(2n + 3)(2n-l)(n + l) 


• no 


® Xo n + 3 Qfim H{n + B) m (®> + 


3n+2 
n + 1 


0nmRnm(0» 


(5-31),. 


+ 


3n + 3 
n-1 


TnmR (n— s) n Rna 


/ 3n + l 
t n + 3 


O^m S (n+ a) It ( 


+ 


3n + 2 
n+1 


finmSnft (®> 


+ 


3n + 3 
n-1 


Tno S^n—2^ B (6, X) ^ • 


It is possible, of course, to write this result in a more convenient form. 
However, we are going to use this result only in order to derive the relationship 
between the disturbing potential on the sphere with radius a from the dislurbing 
potential on the ellipsoid, which is done in the following section. The expression 
(5-31) is very suited for this purpose. 
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6. The Connection Between the PotenMal on the Ellipsoid and the Potential in Space 


fii <3ie previous section we have mapped data from the ellipsoid 
to an auxiliary sphere, solved an integral equation for the desired fimction and 
remapped the solution to the ellipsoid. In this section we will compute ♦■he po- 
tential on the tangent sphere of radius a as a function of boundary values on the 
ellipsoid. 

We can write the potential (of analytical continuation) on the sphere 
with radius a as a series of harmonics 

CO n 

I 

1=C 

complement of the geocentric latitude $ 

R„n(0, X) \ unnormalized spherical harmonics in view of 
Snm(®f the simpler recursion formulae. They have 
the same definition as in (Heiskanen-Morltz, 
1967 ). 

As data on the surface of the ellipsoid we shall consider 6g and later on 
also T and Ag. 

We start the derivation with the representation of the disturbing poten- 
tial outside the ellipsoid by the well-known formula of a surface layer x ’ 


(6-1) T, = y 



( 6 - 2 ) 


T* = 



dE 


together equation fttr the surface layer 


( 6 - 3 ) 


2ttx 


5 

Sn 


J 


X 


k 


dE 


= 6g 
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The definition of the terms, especially of I and can be seen firom the figure, 



direction of the outer ellipsoid normal 

© = 90 - $ 


Similar as in the previous section we must transform the Integral equation for 
the surface layer x mto an equation over the unit sphere, using geocentric lati 
tudes $ as parameters. 

For this purpose, we can use again the ejqjression 


(6-4) 


(—) 

4 - 1 1 

'l -4e'»oos'a' + e<8 (oos@ - oose-f_\ 

\Jk / 

• 2a4 ' 



■taking now the derivative in the jBxed point P’. 

Together with (Molodenskii 1961, page 56) 


(6-5) 


dE 

k 



-f-e'^ sin®$) (1 + -^e’-^ sin^ $') cos $ d$ dX 
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we arrive wif;h the transformation onto fhe sphere at 


(6-6) 


2ttX + * 1 X — + ® 


(T 


t2 




J[- 

a 


— eos^ 0 - cos^'©' 

4 4 


^ (cos© - cos0')^ " I ^ de _ Qg 

ro^ -I -^0 


Now we evaluate the surface layer X in powers of e'‘ 


X = x° + e»®X^ -h e'^x^ 


Substituting this expression in the integral equation (6-6) and equating the coeffi- 
cients of the same power of e'^ we arrive at the following sysiem of equations 


(6-7) 


2ttX° + 


*1 


- 6g 

Ho 


. -L. r 

4tt J 


X 


1 ^ 
Hq 




^J[ 


— cos^0 + cos^©’ 
4 4 


(cos© - cos©')^ 1 „0 

- J X 


do 

Hq 


Again, we are only interested in the first and second equation. With formula 
(5-18) we have also due to the orthogonality relations of spherical harmonics 


(6-8) 


1 

4it 


JXn(®. X) 
a 


dg 

Ad 


1 

2n + 1 


Xn(®’. X’) 
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and with this the solution of the first equation 


(6-9) 


— 7 

4rr Li 

n = f 


2n + 1 
n + 1 


6gn . 


In order to solve the second equation we must expand the right hand side in a 
series of harmonics. As has still to be multiplied by e’® we can use the. 
formulae in spherical approximation 


(6-10) 6g, . T. 

and also 


( 6 - 11 ) 




with 


Then we get 


Tji ~ ^ 1 • 


n = 0 


I 

1 = 1 


2n H- 2 
2n + 1 


1 

4TTa 


{ 


cos^ ©* 
4 



T 


n 


7 V (4m" “ 1) T 

ho ho (2n+3)(2n-l) 
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1 

4 


(2n + 1) . 

(2n + 5) • 


1 i A.. [ 


n =0 B n O 


+ 


Sum Rnn (®» 


^ '..(2p.tiL 
(2n - 3) 



Developing also cos®©' T„ with the help of formula ( 8-5) and comparing the two 
sides of the equation results in 


X 


1 


Sna 


CO n 

{-I I 

n = 0 n=0 


(2n -tl) (4m® - 1) 

(n + 1) (2n + 3) (2n - 1) 


T„„(0, X) 


( 6 - 12 ) 


+ 


CO n 

I I-«[' 


n =0 m =0 


2 (n + 1) 
n + 3 


OnBR(trf3)a(®» ■ + ' 


(2n + 1) 
(n + 1) 


• SnnRnol®* X) + 


- 1 ) _ } 


The solutions (6-9) and (6-12) must be inserted into (6-2). With (Molodenskli 
1962^ page 56) " 


(6-13) 

dE = 

a® (1 - e'® sin® $) cos $ d $ d X 

(6-14) 

li - 

2 • sin((i/2 * /a • rE 

(6-15) 

• 

Te = 

a(l --ce'® sin® $) 

we get 



(6-16) 

dE 

4: ^ (1 - f e'® cos^ 0) cos $ d$ dX 
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and 


(6-17) Ta = J — dE = a J X(1 - -|-e'®cos®0)^cos^ d$ dX 

E GTC ^ 


Evaluating also T into a power series of the small parameter e’*" 


o A. 

Ta = To + e'^ T + e’ -i 


and inserting this expression in (6-17) we obtain 


(6-18) 


T“ = a J X'" 

a 


and 


(6-19) 


P = -aj 


_3_ 

4 


oos^ex” 


^0 5 


We get as an intermediate solution 


( 6 - 20 ) 


T° = a 


n=0 


n + 1 




or written as an integral equation 


(6-21) 


2nT° + i J 

a 


= a 
io 


I ^ 


da 

Hq 


and afterwards 
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T 


1 



n =0 o riO 


(4m^ - 1) 

(n + 1) (2n + 3) (2n - 1) 




( 6 - 22 ) 



I (3n + l,j; A..{^ 

n=0 n=0 ' 


R(n+s)tt(®» 


+ 


En.(®, X) + 

(n + 1) 


(n-l) 


R(n-e)n(®» X) 


} 


Whereas usually the integrations over the unit sphere are made, with the aid of 
the geodetic latitude <f>, the geocentric latitudes $ are used as parameters in 
equation (6-21). - Applying here (Lelgemann, 1970) 


. (6-23) 


cos# = (1 + e'^sin^<p)cos^ 


(6-24) 


(6-25) 


d# = (1 - e’® + 2e'® sin®^) d<p 


1. ^ 1 ^ 1 _ _el! 

JIq 2sin?’/2 2 sin 0/2 L 2 


* (cos® 0 + cos® 6*) 


+ 


(COS0 - cos 9*)® 
4 sin® 0/2 


with 


0 - 90 - ^ 

COS if = cos $ cos f' + sin# sin#’ cos AX 
COS0 = cos<pcos^ + sin^ sin!¥5 cos AX 
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we arrive at the following integral equation system 


T° 


+ 


4n 




+ — i— r = - 

1 

f T° {- 1 -icos 

c 

1 

4n 

J 1 

g 


g 

(6-26) + (COS 0 - cos 0;) 1 

^ ’ 4 sin® 0/2 J 

dg 

%D 

* 2n I®* { 

mm 


cos 


3 0 + (COS e - cos 9')^ \ 

4 sin® ^/2 J 


dg 

As 


fjil3 _ 


with 


da = sin0 d0 dX, 

We have in spherical approximation 

03 


2a 


= I 


n =0 


(2n + 2) T„ 


and therefore 

(6-27) ^ J ^ 

g 


4rr 




g 


® 0» + -|- cos® 0 
~1 - aCOB® 0' 


-I- 1) T„ -1 
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dg 

Ad 


- -Icos^ 0’ + -|- cos® 0 + 


(cos 0 - cos 
4 sin i/j/2 


} 


A series evaluation gives 


(6-28) 


1 _ 1 V (2n + 1) T, , i V ^ (2n H- 1) (4m® - 1) ^ 

“ ® Z. (n + 1) “ ^ L L (n + 1) (2n + 3) (2n - 1) " 

n=o' ^ n=0 u=0 

^ ^ Anm I gnDR(n+S)n(®> + ”^TT** 


n = o m =0 


n - 1 


Adding now (6-22) and (6-28) we get 


= if T + f f • ^ 

n + 1 " A. L (n + l)(2n+ 3)(2n-l) 


n==0 


n = o m=0 


•nffl 


(6-29) 


+ 


■4“ I t { “n +"^ OinR(nva)B(e, X) 


n=.0 n=0 


5n+3„ 5n+7 

n~+~i ^niiRnm(®» X) -t- n— J 'ynnR (o-p)b (®» X)j * 


We sum up the result by 


— 48 — 



cp, n 

T — y' {Ann [ pnmK(&»-g)B (®» Q,{inRaB(®» 

n=o n=0 


+ rnnR(n-g)B (®» ^)] Bnh [ pnaS(t»ta)B(®» '*■ QnjiSnii(6» 


+ rji[jS(ii^^n(0, X)] } 


with 


(5n - l)(n -■ m + 1) (n - m i- 2) 

, (n + 3) (2n + 1) (2n + 3) 

- 6n^ - 8n^ -t 6nm^ >- 6m^ + n + 3 
(n + 1) (2n + 3) (2n - 1) 

(5n H- 7) (n + m) (n + m - 1) 

(n - 1) (2n + 1) (2n - 1) 


In the case of the disturbing potential of the earth we have 

Aoo = Aio == All - Bii ” 0. 

We get a proper expression of the final resxdt from the solution of (6-21) 
together with a simple renumbering of the terms of this series; 

(6-30) Ta = T° + e*® 

CO 

(6-31) (6, X) = a y 6g„(0, X) 

4j n + i 

n — 9 
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(6-32) 


CO “ 

I I 

n=e B=0 


CnaRnti 


(0, A.) + D:„S„a(e,X)l 


Cnm ~ t A^n^ja-pnm -A-nnOttm ■^(n-2)B^nB ] 


” EB^n-EjnPnn BnmQnn + B^-a^arnnl 


With 


(5n - 11) (n - m - 1) (n - m) 

4 (n + 1) {'211 - 3) (2n - 1) 

- 6n^ - 8n^ )• 6nm^ - 6m^ + n h 3 
4(n I- l)(2n-l)(2m-3) 

SE (5n + 17) (n + m + 2) (n 4- m -i- 1) 

4(n + l)(2n -1- 5)(2n + 3) ' 


It may be pointed out that the gravity disturbances belong to points Pe {(p, X) on 
the ellipsoid and the disturbing potential to points P^ (5, X) on the sphere. 
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Based ob this result, we will derive some other useful formulae. We begin 
with a direct relation between the gravity disturbances at the ellipsoid 6gE 
and the gravity disturbances ^ at the tangent sphere 


(6-33) 


6gk = - 6g3 . 


We use the abbreviation (Sgs) because the abbreviation (Sgs) is already 
reserved for the gravity anomalies (see formula (4-6) ). The potential T» 
is given on a sphere. So we have the exact series relation 


(6-34) I (n + l)(Ta)«. 

n=0 

hiserting (6-30) into (6-34) we get as the desired result 


(6-35) 


2 ® n 

6gj = 6gE + y (n + 1) V [ CnaRno(®» "*■ RnaSnii(6, X) ]. 


n = 0 


D=0 


or 


y (n + 1) y [ C^BRnn(t^» X) + p^S{im(y, X) ] , 

n=: 0 m= 0 


Next, we will derive a direct relation between the disturbing potential T, at the 
tangent sphere and the disturbing potential Te at the ellipsoid. This derivation 
should be performed in three steps. First we compute from T* with the 
help of (6-20). Second, we compute the gravity disturbances at the ellipsoid with 
(6-33), using for more convenience the intermediate .relation (6-29), Third, we 
compute the disturbing potential Te at the ellipsoid with formulae (5-19) and 
(5-31) of the previous section. We arrive with the simple result 


(6-36) 


Te - Ta + 6t3 
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with 


oo n 

6ta = — ~ — ^ { Ann (®» ■*■ 

n = S n=0 


+ ynnll{n=a)!» (®» 1 ^no E 0oaS(nh5)m(®» ^) 


+ ftia Rnn ( ^) ■*■ T no (n-a)o { 1 } ’ 


or with the aid of ( 8-5) 


(6-37) 6ts = -|^ cos® e T (6, X ) . 

Because Stg is a term of second order. we can neglect the difference between 
6 and ©, writing also 


Sts 



cos® 0 T (0, X) . 


Finally, we shall derive the relation between the gravity anomalies at the 
tangent sphere and the gravity anomalies Agg at the ellipsoid, s<^arting from 
the known relation at the ellipsoid (Lelgemann, 1970) 


(6-38) AgE = 6g + ^(1+e'®- e’®cos®6)TE. 

/ on 3. 

We have by definition 

(6-39) Agi = AgE - 6gs 
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and 


2 

(6-40) Ag, = - — Ta. 

After some manipulation we get from these formula 

6g3 =''6g3 - — (1 +e’2 - e'^cos^6)TE + ™Ta 
a S' 

or as a final result 

(6-41) 6 g 3 =■ figs - (2 - -r cos^ 0) T ( 6 , X) 

ci 

It is very easy to derive (4-16), (4-17), (4-18) and (4-19) from (6-33), (6-35) 
and (6-41). 
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7. Summary 


The definition of the geoid is discussed in the first part of the present 
study. Of course, the geoid is a certain equipotential surface of the earth gra- 
vity field. The discussion is only concerned with a specification of this equipo- 
tential surface. The final proposition for a definition is based on the geodetic 
height systems by an inclusion of given information from oceanographic science 
about sea surface topography. 

The realisation of the geoid should be made in two steps. First, ehe 
existing height systems are corrected by information about sea surface topo- 
graphy, An equipotential surfe.ce is choosen in such a manner that the sum of 
the squares of the deviations to the corrected main height systems is a mini- . 
mum. This particular equipotential surface is called the geoid. 

A least squares procedure is derived for a realisation of this definition. 

The precise cartesian coordinates of geophysical stations and in addition the 
quasi geoidal undulations and the normal heights in these points are needed as 
data. The unknown corrections to the various height datums influence, of course 
the gravity anomalies and consequently ihe geoidal undulations. This influence 
is regarded in the least squares procedure. 

In the second part of the study a solution of Stokes* problem with an ac- 
curacy of better5than ±10 cm is suggested, in which the results of satellite geodesy 
can be' included in a rather simple way. The method is based on the possibility that 
the potential outside the earth surface can be approximated by the potential of ana- 
lytical continuation inside the earth with any accuracy. Gravity anomalies from the 
earth’s surfece are reduced by three successive corrections to gravity anomalies at 
the sphere with radius a. At this sphere they can be combined with given potential 
coefficients from satellite geodesy. Then the potential at the surface of the earth 
is computed with the aid of additional three corrections. 

Of special interest could be the very simple formula which connects 
the potential at the ellipsoid with the potential at the sphere with radius a. In 
the computation of potential coefficients from altimeter data ^he influence of 
ellipticity can easily be taken into account with this formula. 
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8. Appendix 


We shall derive in the appendix an evaluation of the function 

(cos 9 - cos 6')^ 

8 sin^ 0/2 

6. , . . polar distance 

. spherical distance between two points 
P' and P 


in a series of spherical' harmonics, starting with the derivation *of two formulae 
which we shall need in the final evaluation. 

As the immediate result of the kernel of Poisson’s integral formula 
(Heiskanen-Moritz, 1967) 

(2n+l) (-^Jp„(cos^[)) 

n=n ^ 


we get for r = R 

00 

(8-2) ^ (2n + l)P„(cos0) = 0. 

n =0 

Now, we will try to deriveia series ev^uation for 


1 ^ q 

8sixflj)/2 


In this case^ we do not get a well defined expression due to the relation of (8-2). 
The evaluation can be started with the well known formula 
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jj [ (r - + 4 r R sin^ !^/2 r 




The first' derivative gives 


aqA) ^ 

. 3r 


CO 

-I. 


(n + 1) 


R” 


n = 0 




PqcoS!/)) 


t(r - R) + 2Rsin® ip/2) 


The second derivative gives the result 


= y (n+l)(n + 2) Pn(cos0) 

Sr n^O 


^ + 3 [ (r-R) + 2Rsin^0/2] 


3 ill /o 1 2 


On the sphere with the radius R we have 


£o = 2R sin 0/2 


For points on this sphere we can write 


^ y (a + 1) (n + 2) Pn(cos0) 

ri“0 


8R^sin^0/2 . ^ 8 R® sin 0/2 
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A further evaliaation gives 


CO 

SB^W ‘ ^ + - 


n=:0 


4 2 sin 0/2 


or 


CO 

Sslhir = I [(n + l)(n + 2)--|l] P,(oos#) 


n = 0 


With the identity 


(n + l)(n+2)- -f 


-i- (2n + 1) (2n + 5) 


we get 


1 

8 sin^ 0/2 


00 

^ (2n + 1) ( 2n + 5) Pn(cos0). 

n=; 0 


Because of (8-2) we can see that the evaluation of l/l^ inio a series of 
spherical harmonies gives not a well defined result but an expression of the 
form 


(8-3) 


1 

8 sin^ 0/2 


I P.icoB 

n zz O 


where a is any given number. On the other hand we shall see that in the 
series expression for formula (8-1) the term a drops out. 
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Let us si-art with the following expression for a fixed point P' (Heiskanen- 
Moritz, 1967). 


^ ^ X E a„=(0',X’)Rn„{e, X) 

^ n 0 111 r= 0 

+ b„„(0', X') S„„(0, X)]. 

Then we can estimate the coefficients ann by 

• a = - ^ (n - m)_i_ [ { wcos 6 - cos 0')^. 

2tt (n + m) ! J ^ ^ ^ 


n = 0 “ 


da. 


Using twice the recursion formula for Legendre functions 


cos 6 Pnj,(cos0) 


n - m + 1 
2n + 1 


R(tt+l) m (COS0) + ■ ” ^ Pt._,l„{cos0) 

2n H- 1 


and multiplying the resuP with cos (mX) we get the following two expressions 
(S-4) cos 0 R„,{0, X) = R(n+i)»(e. X) + ~ R(«-^ »(0. X) 

and 

(S^S) COS^0R„5(0, X) = C^mR^n+s)D(®» X) + 0nitiRniii(®» X) + Xna R(it-a)m(®* X) 
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with 


(n - m -f- 1) (n - m i- 2) 
■(2n + 1) (2h -I- 3) 

2n^ - 2m^ + 2n - 1 
(2n + 3) (2n - 1) 

(n + ro) (n + m - 1) 

■ ( 2n + 1) (2n - 1) 


With the help of this recursion, formula we get, due to the orthogonality relation® 
of spherical harmonics, after the integration 

a„„(0«, x>) = [ -(2n + 1) (n + a -i- 2) : 

(n -t- m) ! ' ' 

- (2n + 1) (n I- a) X’) 

-.(2n + 1) (n + a - '2) ynoR X') 

+ 2 cos 9* (n " m + 1) (n + a ■!• 1) R(n+i)ii'(®'» X’) 

+ 2 cos 9' (n + m) (n + a ~ 1) R(»^) m (9'» X') 

- cos®6'(2n -tl)(n +a)R„o(e', X*) : 


Using the recursion formula once more the term a drops out and we are left 
with the result 


= - 2 ~ ! r (n - m -f 1) (n + m + 1) 

(n + m) ! L . (2n i 3) 


(n -t m) (n » m) 
(2n - 1) 




• Rn»(0'. X') 
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or 


a. 


4x11.^ 1 

(n + m) 1 (2n-f3) (2n - 1) 


2 


Analogoxis evaluations give 


^ = 2 <° - g)-'- 

(n + m) I 


4m^ - 1 


(2n + 3) (2n - 1) 


S„b(0’, X*). 


As a final result, we get the desired expression 


(8-6) 


(cos 9 - cos 0')^ 
8 sin^ ip/2 


CQ 

I {- 


n = 0 


(2n + 3) (2n - 1) 


P„( cos 0) P„(cos 0') 


+ 2 


V M 

A fn 


" m) ! 


4m^ - 1 


(n + m) ! (2n + 3) (2n - 1) 

+ S„40', X>) S„„(0, X) ] }. 


[Rn„(0\X’)R„„(0, X) 
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